Spinor Dynamics in an Antiferromagnetic Spin-1 Condensate 
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We observe coherent spin oscillations in an antiferromagnetic spin-1 Bose- Einstein condensate 
of sodium. The variation of the spin oscillations with magnetic field shows a clear signature of 
nonlinearity, in agreement with theory, which also predicts anharmonic oscillations near a critical 
magnetic field. Measurements of the magnetic phase diagram agree with predictions made in the 
approximation of a single spatial mode. The oscillation period yields the best measurement to date 
of the sodium spin-dependent interaction coefficient, determining that the difference between the 
sodium spin-dependent s-wave scattering lengths a/=2 — a/=o is 2.47 ± 0.27 Bohr radii. 

PACS numbers: 03.75.Mn, 32.80.Cy, 32.80.Pj 



Atomic collisions are essential to the formation of Bose- 
Einstein condensates (BEC), redistributing energy dur- 
ing evaporative cooling. Collisions can be coherent and 
reversible, leading to diverse phenomena such as super- 
fluidity [lj and reversible formation of molecules [2| in 
BECs with a single internal state. When internal de- 
grees of freedom are included (as in spinor condensates) , 
coherent collisions lead to rich dynamics [H, 3 in which 
the population oscillates between different Zeeman sub- 
levels. We present the first observation of coherent spin 
oscillations in a spin-1 condensate with antiferromagnetic 
interactions (in which the interaction energy of colliding 
spin-aligned atoms is higher than that of spin-antialigncd 
atoms.) 

Spinor condensates have been a fertile area for the- 
oretical studies of dynamics ground state 
structures d, [l(| , and domain formation 11 1 . Extensive 
experiments on the ferromagnetic F—l hyperfine ground 
state of 87 Rb have demonstrated spin oscillations and co- 
herent control of spinor dynamics [H, EH • Observation of 
domain formation in 23 Na demonstrated the antiferro- 
magnetic nature of the F—l ground state [l]| and de- 
tected tunneling across spin domains 1J| ; no spin oscilla- 
tions have been reported in sodium BEC until now. The 
F—2 state of 87 Rb is thought to be antiferromagnetic, but 



a cyclic phase is possible [la, LL6j • Experiments on this 



state have demonstrated that the amplitude and period 
of spin oscillations can be controlled magnetically [4j . 

At low magnetic fields, spin interactions dominate the 
dynamics. The different sign of the spin dependent in- 
teraction causes the antiferromagnetic F=l case to differ 
from the ferromagnetic one both in the structure of the 
ground-state magnetic phase diagram and in the spinor 
dynamics. Both cases can exhibit a regime of slow, an- 
harmonic spin oscillations; however, this behavior is pre- 
dicted over a wide range of initial conditions only in the 
antiferromagnetic case [8] . The spin interaction energies 
in sodium are more than an order of magnitude larger 
than in 87 Rb F=l for a given condensate density [3|, 
facilitating studies of spinor dynamics. 



The dynamics of the spin-1 system are much simpler 
than the spin-2 case [1, HH, [l|| , having a well-developed 
analytic solution Q. This solution predicts a divergence 
in the oscillation period (not to be confused with the 
amplitude peak observed in 87 Rb F=2 [4j oscillations). 

This Letter reports the first measurement of the 
ground state magnetic phase diagram of a spinor con- 
densate, and the first experimental study of coherent 
spinor dynamics in an antiferromagnetic spin-1 conden- 
sate. Both show good agreement with the single-spatial- 
mode theory [lfj. To study the dynamics, we displace 
the spinor from its ground state, observing the resulting 
oscillations of the Zeeman populations as a function of 
applied magnetic field B. At low field the oscillation pe- 
riod is constant, at high field it decreases rapidly, and at a 
critical field it displays a resonance-like feature, all as pre- 
dicted by theory 0. These measurements have allowed 
us to improve by a factor of three the determination of 
the sodium F = 1 spin-dependent interaction strength, 
which is proportional to the difference a/=2 — <2/=o in 
the spin-dependent scattering lengths. 

The state of the condensate in the single-mode ap- 
proximation (SMA) is written as the product (f>(r)( of a 
spin-independent spatial wavefunction (f>(r) and a spinor 
C = {JpZe ld ~ ,^Me ie °,^/p^e ie +). We use p_, p , and 



p + (9-, 6*o, and 6+) to denote fractional populations 
(phases) of the Zeeman sublevels mp = —1, 0, and 1, 
so that Pi = l- The spinor's ground state and its non- 
linear dynamics may be derived from the spin-dependent 
part of the Hamiltonian in the single-mode and mean- 
field approximations, subject to the constraints that to- 
tal atom number N and magnetization m = p+ — p~ are 
conserved 0. The "classical" spinor Hamiltonian E is a 
function of only two canonical variables: the fractional 
population po and the relative phase 8 = 8 + + 8- — 28q. 
It is given by 

E = 5(1- po) + cp ((l-p ) + \J (l-pof-m 2 cosd^j, 

(1) 

where 5 = h x (2.77 x 10 10 Rz/T 2 )B 2 is the quadratic 
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Zeeman shift [8| with h the Planck constant. (The linear 
Zeeman shift has no effect on the dynamics.) The spin- 
dependent interaction energy is c — c 2 (n) , where (n) is 
the mean particle density of the condensate and 



(>2 



3M 



(a/ =2 - a/=o) 



(2) 



is the spin-dependent interaction coefficient [8|, iTjj . Here 
M is the atomic mass, at— 2 and a/ = o are the s-wave 
scattering lengths for a colliding pair of atoms of total 
spin / = 2 and / = 0, respectively; Bose symmetry en- 
sures there are no s-wave collisions with total spin of 
1. If C2 is positive (negative), the system is antiferro- 
magnetic (ferromagnetic). The spinor ground state and 
spinor dynamics are determined by Eq. ([T]). 

The ap paratus is similar to that described previ- 
We produce a BEC of 10 5 
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Na atoms in 



1 state, with an unobservably small thermal frac- 
tion, in a crossed-beam 1070 nm optical dipole trap. The 
trap beams lie in the horizontal xy plane, so that the trap 
curvature is nearly twice as large along the vertical z axis 
as in the xy plane. By applying a small magnetic field 
gradient with the MOT coils (less than lOmT/m) during 
the 9 s of forced evaporation, we fully polarize the BEC: 
all atoms are in mp = +1. Conservation of spin angular 
momentum ensures that the magnetization remains con- 
stant once evaporation has ceased; a state with p + = 1 
persists for the lifetime of the condensate, about 14 s. 

We then turn off the gradient field and adiabatically 
apply a bias field B of 4 to 51 /iT along x, leaving the 
BEC in the p + = 1 state. To prepare an initial state, 
we apply an rf field resonant with the linear Zeeman 
splitting; typically the frequency is tens to hundreds of 
kilohertz. Rabi flopping in the three-level system is ob- 
served 19], and controlling the amplitude and duration 
of the pulse can produce any desired magnetization m, 
which also determines the population po. The flopping 
time is less than 50 /is, much shorter than the character- 
istic times for spin evolution governed by Eq. |T]). Using 
this Zeeman transition avoids populating the F=2 state, 
thus avoiding inelastic losses, which are much greater for 
23 Na than for 87 Rb. 

We measure the populations pi of atoms in the three 
Zeeman sub-levels by Stern-Gerlach separation and ab- 
sorption imaging [20j . The Stern-Gerlach gradient is par- 
allel to the bias field £?, while the imaging beam propa- 
gates in the i direction. The phase 9 is not measured. 

To measure the ground state population distribution 
as a function of magnetization and magnetic field, we 
first set the magnetization using the rf pulse. We then 
ramp the field to a desired final value over 1 s, wait 3 s 
for equilibration, and measure the populations as above. 

Figure QJb) displays the measured ground-state mag- 
netic phase diagram. The theoretical prediction in 
Fig. [lja) is the population po that minimizes the energy, 




FIG. 1: a) Theoretical prediction of the ground-state frac- 
tional population po as a function of magnetization m and 
applied magnetic field B, assuming a spin-dependent interac- 
tion energy c — /ix20.5 Hz. The thick line lying in the po — 
plane indicates the boundary between the po — and the 
po > regions, b) Experimental measurement. The surface 
plot is produced by interpolation of data points. 



Eq. ([T]) . Such minima always occur at 9 = tt for antifer- 
romagnetic interactions. The measurements agree well 
with the prediction, which is made for spin interaction 
energy c — hx 20.5 Hz (determined by spin dynamics as 
described below). 

The first term of Eq. |T]) depends on the external mag- 
netic field and tends to maximize the equilibrium po 
population. The second, spin dependent, term has the 
same sign as C2 and in the antiferromagnetic case tends 
to minimize the equilibrium po population. The phase 
transition indicated by the thick line in Fig. [IJl arises 
at the point where these opposing tendencies cancel for 
po = 0. Along the transition contour, po rapidly falls 
to zero. By contrast, the ferromagnetic phase diagram 
has po = only at m = 1. In the region B < 15 pT 
and m > 0.6, there should be virtually no population in 
mp = for antiferromagnetic interactions, and popula- 
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tions up to po — 0.34 for ferromagnetic interactions (as- 
suming the same magnitude of c). For our equilibrium 
data, the reduced % 2 with respect to the antiferromag- 
netic (ferromagnetic) prediction in this region is 2 (20). 
This demonstrates that sodium F = 1 spin interactions 
are antiferromagnetic, as previously shown by the misci- 
bility of spin domains formed in a quasi-one-dimensional 
trap [13| . 

Across most of the phase diagram, the scatter in the 
population is consistent with measured shot-to-shot vari- 
ation in atom number. This variation is 20%, implying 
an 8% variation in the mean condensate density accord- 
ing to Thomas-Fermi theory. The variance of results is 
not due to the magnetic field (calibrated to a precision of 
0.2 /LtT), nor to residual field variations across the BEC 
(less than 250 pT). Uncertainties in setting the magne- 
tization are obviated, as the magnetization is measured 
for each point as the difference in fractional populations 
m = p+ — p- . Discrepancies between theory and exper- 
iment at low magnetic fields may be attributed to the 
field dependence of the equilibration time. We observe 
equilibration times (see below) ranging from 200 ms at 
high fields to several seconds at low fields, by which time 
atom loss is substantial. 

If the spinor is driven away from equilibrium, the full 
coherent dynamics of the spinor system Eq.([T]) are re- 
vealed. We initiate the spinor dynamics with the rf tran- 
sition described above, but now look at the evolution over 
millisecond timescales. 

The spinor dynamics are described by the Hamilton 
equations for Eq. (p]) Q: 

2 dE ■ 2 BE 

The system is closely related to the double- well "bosonic 
Josephson junction" (BJJ) [2l|, [22| and exhibits a regime 
of small, harmonic oscillations and, near a critical field 
B c , is predicted to display large, anharmonic oscilla- 
tions. At B c the period diverges (where 8(B C ) = c[(l — 
po) + \/(l — Po) 2 — m2 cos 9} , with po and 9 taken at 
t = 0) 8]. The critical value corresponds to a transi- 
tion from periodic-phase solutions of Eq. ((3|) to running- 
phase solutions. At the critical value it is predicted that 
the population is trapped in a spin state with pa =0. This 
phenomenon is related to the macroscopic quantum self- 
trapping that has been observed in the BJJ [22j]. How- 
ever, very small fluctuations in field or density will drive 
Po away from 0. Observing a ten-fold increase in the pe- 
riod above its zero-field value would require a technically 
challenging magnetic field stability of better than 100 fT. 

Figure [5] plots the period and amplitude of oscillation 
as a function of magnetic field. An example of the os- 
cillating populations is shown in the inset. The spinor 



condensate is prepared with initial po — 0.50 ±0.01 1 and 
rn = 0.00 ± 0.02, and a plot of p versus time is taken 
at each field value. Qualitatively, the period is nearly 
independent of magnetic field at low fields, with a small 
peak at a critical value B c = 28 pT, followed by a steep 
decline in period. The amplitude likewise shows a max- 
imum at B c . Oscillations are visible over durations of 
40 ms to 300 ms. Beyond these times, the amplitude of 
the shot-to-shot fluctuations in po is roughly equal to 
the harmonic amplitude. This indicates dephasing due 
to shot-to-shot variation in oscillation frequency, proba- 
bly associated with the variations in magnetic field and 
condensate density, rather than any fundamental damp- 
ing process. At even longer times, we observe damping 
and equilibration to a new constant po\ the damping time 
varies with magnetic field from 200 ms to 5 s. 

For the theoretical prediction in Fig. [2) the initial value 
of po and m are obtained experimentally. We treat only 
c and 9(t — 0) as free parameters; c is also predicted 
by prior determinations of C2 and our knowledge of the 
condensate density. The initial relative phase is not the 
equilibrium value 9 = tt, due to our rf preparation. For a 
three-level system driven in resonance with both transi- 
tions, the relative phase is 9 = at all times during the rf 
transition, as we derive from Ref. [l9| . Small deviations 
from initial 9 = could be caused by an unequal splitting 
between the levels, from e.g., the quadratic Zeeman shift. 

The best fit to the data in Fig. [2^ and b is obtained 
by using c=hx(21 ± 2) Hz and 0(t = O) = 0.5 ± 0.3 (with 
no other free parameters). Away from the critical field 
B c , agreement with theory is good. The fitted value of 
c implies that B c is 27 /iT, in reasonable agreement with 
the apparent peak observed at 28 pT. Our ability to ob- 
serve strong variations in period near B c is limited by 
density fluctuations (8%) and magnetic field fluctuations 
(0.2 /iT). Near B c , typically only one cycle is visible be- 
fore dephasing is complete. Such rapid dephasing can, 
itself, be taken as evidence of a strongly independent 
period, as expected near the critical field. 

To include the known fluctuations in density and mag- 
netic field in our model, we perform a Monte Carlo sim- 
ulation of the expected signal, based on measured, nor- 
mally distributed shot-to-shot variations in values of c, 
S, m and po(t = 0). At each value of B in Fig. [51 we 
generate 80 simulated time traces, with each point in the 
time trace determined from Eq. [3] We fit the simulated 
traces using sine waves and record the mean and stan- 
dard deviation of the amplitude and period of the fits. 
The results (shaded regions in Fig. [2]) show a less sharp 
peak in the period. The smoothing of the peak at B c is 
consistent with our data. 



All uncertainties in this paper are one standard deviation com- 
bined statistical and systematic uncertainties 
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FIG. 2: Period (a) and amplitude (b) of spin oscillations as a 
function of applied magnetic field, following a sudden change 
in spin state. The solid lines are theoretical predictions from 
solving Eq. ([3]). The theoretical prediction of the period goes 
to infinity at about 27 pT. The shaded regions are ±1 stan- 
dard deviation about the mean values predicted by the Monte 
Carlo simulation. Inset: Fractional Zeeman population (solid 
dots) and magnetization (open circles) as a function of time 
after the spinor condensate is driven to po = 0.5, m = 0. 
B = 6.1 fiT. The solid line is a sinusoidal fit. 



It is clear in Fig.[2]that the oscillation period is insensi- 
tive to the magnetic field at low values of the field. In this 
regime, the period is sensitive only to the spin interaction 
C2 and the density of the condensate (n) . Measuring this 
period allows us to determine the difference in scattering 
lengths a/ = 2 — a/=o- The trace inset in Fig. [5] was taken 
in this regime, at a magnetic field of B = 6.1 p;T, and 
shows harmonic oscillations with period 24.6 ± 0.3 ms. 
Here the predicted period dependence on magnetic field, 
is indeed weak and the oscillations dephase 
only slightly over the duration shown. Using this mea- 
surement of the period (in which much more data was 
taken than for each point making up Fig. [5] (a) and (b)), 
and including uncertainties in initial 8, po, and m, we 
obtain the spin interaction energy c— h x (20.5 ± 1.3) Hz. 

Finding a/ = 2 — a f=o requires a careful measurement 
of the condensate density. We take absorption images 
with various expansion times to find the mean field en- 
ergy. The images yield the column density in the xy 
plane, and the distribution in the z direction can be 
inferred from our trap beam geometry. We find that 
the mean density of the condensate under the conditions 



of the inset to Fig. [2] is (n) = 8.6 ± 0.9 x 10 13 cm~ 3 . 
From this we calculate a/ = 2 — a f=o = (2.47 ± 0.27)ao, 
where an = 52.9 pm is the Bohr radius. This is consis- 
tent with a previous measurement, from spin domain 
structure, of a/=2 — a f=o = (3-5 ± 1.5)an [HI and is 
smaller than the difference between scattering lengths 
determined from molecular levels, a/ = 2 = (55.1 ± 1.6)an 
and a/ = o = (50.0± 1.6)ao [23J]. A multichannel quantum 
defect theory calculation gives a/ = 2 — a/=o = 5.7an Q. 

Finally, we consider the validity of the spatial single- 
mode approximation. The SMA was clearly violated in 
previous work on 23 Na [l3| and 87 Rb Q F—l spinor 
condensates where spatial domains formed. Spatial de- 
grees of freedom decouple from spinor dynamics when the 
spin healing length £ s = 2irh/ \/2m\c2\n is larger than the 
condensate. From our density measurements we find typ- 
ical Thomas-Fermi radii of (9.4, 6.7, 5.7) /im. The spin 
healing length, based on our measurements of c, is typi- 
cally t; s = 17 /im. We therefore operate within the range 
of validity of the SMA. Furthermore, Stern-Gerlach ab- 
sorption images show three components with identical 
spatial distributions after ballistic expansion, indicating 
that domain formation does not occur. 

In conclusion, we have studied both the ground state 
and the spinor dynamics of a sodium F=l spinor conden- 
sate. Both agree well with theoretical predictions in the 
SMA. By measuring the spin oscillation frequency at low 
magnetic field, we have determined the difference in spin- 
dependent scattering lengths. The observed peak in oscil- 
lation period as a function of magnetic field demonstrates 
that the spinor dynamics are fundamentally nonlinear. 
It also suggests the existence of the predicted regime of 
highly anharmonic spin oscillations at the center of this 
peak, which should be experimentally accessible with suf- 
ficient control of condensate density and magnetic field. 
Observation of anharmonic oscillations, as well as popu- 
lation trapping and spin squeezing effects, could be aided 
by a minimally destructive measurement of Zeeman pop- 
ulations 25[ to reduce the effects of magnetic field drifts 
and shot-to-shot density variations. 
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